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SURFACE-REFRACTIVITY MEASUREMENTS 
AT NASA SPACECRAFT TRACKING SITES 


by 

P. E. Schmid 

Goddard Space Flight Center 


INTRODUCTION 

Tracking data bias caused by the refractivity of the troposphere (and the remainder of the un- 
ionized atmosphere) is elevation-angle dependent: At 5° elevation, the bias may be up to 20 m in 
range, 50 cm/s in range rate, and several milliradians in angle (Reference 1); at the zenith, the tropo- 
spheric range bias decreases to 2 or 3 m, and the range rate and angle tropospheric biases tend to zero. 
The ionized portion of the upper atmosphere, the ionosphere, can also be a significant bias source, 
especially at spacecraft tracking frequencies below 2000 MHz. However, ionospheric effects are not 
discussed in this paper. 

The frequency-independent tropospheric refraction effect on angle of arrival data is highly cor- 
related with surface refractivity at the reception point of the radio wave. Range and range rate data 
tropospheric biases are due to integrated propagation delay effects; the greatest source of variability is 
the “wet term” contribution in the near-Earth refractivity. Refractivity N is simply a convenient 
means for specifying index of refraction n; the relation between the two is given by 

n = 1 + N( 1 O' 6 ) . 

In this report, the measured time-dependent variation in surface refractivity is presented for 1 2 
Apollo sites during 2 months in 1967 and for 15 sites during 8 consecutive months in 1968. This 
information is useful in estimating the absolute tracking accuracy attainable if, as generally is the case, 
average values of surface refractivity are used in the reduction of tracking data from any given site. 
The original purpose of this test was to ensure that the estimated monthly surface-refractivity values 
used by NASA’s Manned Spacecraft Center (MSC) were adequate for the Apollo program. It was seen 
that Apollo tropospheric corrections were consistent with the spacecraft tracking accuracy required 
for the mission. All Apollo stations were assigned monthly average values of refractivity which in no 
case differed by more than 1 3 percent from the measured monthly averages. 
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However, for tracking experiments such as very-long-baseline interferometry (References 2 and 3), 
the tropospheric variations throughout the day or month cannot be neglected. This is also true for the 
stringent tracking accuracy associated with geodesy and Earth-physics studies (Reference 4), such as 
those now being undertaken at Goddard Space Flight Center (GSFC). 

The best correction for tropospheric bias would of course be obtained from a measured profile 
along the ray path between tracking site and spacecraft. However, such information is not usually 
available; instead, some type of mathematical model is normally assumed. For purposes of radio 
refraction and time-delay correction, an exponential troposphere such as that presented by the 
National Bureau of Standards offers good agreement with measured profiles for most tracking appli- 
cations (References 1 and 5). The scale height for such a model at any given site can be determined 
empirically as a function of surface refractivity (Reference 5) or derived from knowledge of the sur- 
face refractivity and the vertical integrated profile, the dry component of which is only a function of 
surface pressure. 

It should be noted that for short-baseline interferometers, such as in the NASA Minitrack system, 
there is an inherent first-order correction for tropospheric biasing effects. This inherent correction 
results whenever one uses the free-space value of the wavelength of the arriving radio signal to calcu- 
late the short-baseline interferometer angle (Reference 6). 

There are, of course, other tropospheric errors which may predominate, especially at elevation 
angles below about 10°. Horizontal gradients are not included in the usual model, in which a spheri- 
cally stratified (i.e., horizontally homogeneous) troposphere is assumed. Also, at low elevation angles, 
tropospheric turbulence can introduce noise into the tracking measurements. However, these noise- 
like errors tend to be averaged out by the data-smoothing and orbit-computation processes and are 
therefore less serious than steady-state or biasing errors. It should be noted that the noise considered 
in the foregoing sense is not a radio noise, such as that attributed to thermal emission by the atmos- 
phere, but rather a phase modulation of radio signals traversing regions of time-varying permittivity 
and hence time-varying radio index of refraction. 

CALCULATION OF SURFACE REFRACTIVITY 

The relationship between surface refractivity and total tropospheric biasing effects on tracking 
data is well established (Reference 7). This report is concerned only with presenting the means for 
calculating surface refractivity based on surface meteorological data and showing how this value of 
refractivity varies from site to site throughout the year. The relationship between meteorological data 
and refractivity is also well known. The relationship reported by Smith and Weintraub in 1953 
(Reference 8) is still in widespread use by radio scientists throughout the world (Reference 9). Ref- 
erence 9 presents an excellent review of the physics relating the radio dielectric constant (and, there- 
fore, index of refraction) to the parameters of temperature, total air pressure, and relative humidity; 
the derivations will not be repeated here. 

Scope of Test 

Meteorological data were collected from the 1 5 Manned Space Flight Network (MSFN) stations 
listed in Table 1. Station personnel measured the air temperature, wet-bulb temperature (using various 
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Table 1 — MSFN sites that participated in the surface-refractivity test. 


Station 

Station Name 

Latitude 

(deg) 

Longitude 

(deg) 

Elevation Above 
Mean Sea Level 
(m) 

MIL 

Merritt Island, Florida 

28.50 8 272 

-80.693 417 

9.2 

BDA 

Bermuda 

32.351 286 

-64.658 181 

22.6 

GBM 

Grand Bahama 

26.632 858 

-78.237 664 

11.4 

ANG 

Antigua Island 

17.016 692 

-61.752 689 

34.4 

CRO 

Carnarvon, Australia 

-24.907 592 

113.724 247 

44.5 

HAW 

Hawaii 

22.124 897 

-159.664 989 

1150.9 

GYM 

Guay mas, Mexico 

27.963 206 

-110.720 850 

23.9 

TEX 

Corpus Christi, Texas 

27.653 750 

-97.378 469 

12.3 

GWM 

Guam 

13.309 244 

144.734 414 

92.1 

GDS 

Goldstone, California 

35.341 694 

-116.873 289 

973 

ACN 

Ascension Island 

-7.955 056 

-14.327 578 

544.2 

HSK (CNB)* 

Canberra, Australia 

-35.597 222 

148.979 167 

1130.8 

MAD 

Madrid, Spain 

40.455 358 

-4.167 394 

785.1 

CYi 

Grand Canary Island 

27.764 536 

-15.634 536 

160.4 

TAN 

Tananarive, Madagascar 

-19.018 055 

47.304 444 

1385.2 


*Call letters changed during test to HSK (Honeysuckle Creek). 


| types of psychrometers), and air pressure at an average time interval of 6 hours. The test 1 was con- 
ducted throughout August and September 1967 and then again on a continuous basis from February 
( 1 968 to September 1 968. The raw data were mailed to GSFC from each station and entered onto 

computer cards which were then used as input to a computer program. This program performed the 
| appropriate conversion of units, refractivity calculations, and plotting. The results were excellent: 
Based on the reasonable assumption that readings were made to within ±0.25 K in temperature, the 
| absolute error in the surface refractivity N is less than 3 A'-units (Reference 9). The corresponding 
relative changes (i.e., resolution) are of the order of ±0.3 Akmits. Thus, the accuracy of the individual 
refractivity values calculated in this report is of the order of 1 percent (since N $ is of the order of 300). 
The corresponding resolution is 0. 1 percent. 


^P. E. Schmid, “MSFN Site Meteorological Data,” memorandum to R. V. Capo, Manned Flight Operations Branch, GSFC, June 19, 1967. 


3 







Mathematical Formulation of Refractivity 


All data mailed from the tracking sites included the time associated with each measurement. 

Since several different types of psychrometers, thermometers, and barometers were used by the various 
stations, some data conversion was required. Where necessary, the following relationships were used; 


p 

r (mb) 

T CC) 

T (K) 


P r . , X 33.863 95 , 

(inches mercury) 5 


(7>f) * : 32) x 9 , 


T Cc) +273.16 . 


The tropospheric surface index of refraction n s , which is greater than unity, can be written as 


n s = 1 +N s ( 1(T 6 ). 


( 1 ) 


The value of N ranges from a minimum of 200 to a maximum of 450, depending on geographical 
location, time of year, weather conditions, and so on (Reference 10) since N s is a function of atmos- 
pheric temperature, pressure, and water-vapor content (Reference 9). That is, 


77.6 

T 



48 1 Qg\ 

T ) ' 


( 2 ) 


where T is the air, or dry-bulb, temperature (Kelvin);/’ is the total air pressure (in millibars); and e is 
the partial water-vapor pressure (in millibars). In turn, e is related to psychrometric measurements by 
(Reference 9) 

<? = <?' -0.000 61 {T - T)P , (3) 


where e' is the saturation vapor pressure (in millibars) at T' and T' is the wet-bulb temperature 
(Kelvin). The value of e' can be obtained directly from appropriate psychrometric charts (Relerence 
1 1). However, in order to facilitate computer computation, an empirical formulation was used, which 
resulted in negligible error relative to the uncertainty in the meteorological field measurements. The 
expression used to calculate the saturated vapor pressure in millibars is given by (Reference 12) 

e' s =6.11 X io flf / ( * +f) , (4) 

where a = 7.5, t is the wet-bulb temperature (Celsius), and b = 237.3. Note that 

T' = t+ 273.16 (Kelvin). 

Thus, the value of N , which can be directly related to radio-wave tropospheric refraction and time 
delay, is determined by three basic measurements: (1) dry-bulb temperature, (2) wet-bulb tempera- 
ture, and (3) total atmospheric pressure. 
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Once all the data for a particular month from a given station were processed, the monthly mean 
and standard deviation (SD) were computed, with 


mean = N = 


| m - 

where m is the number of observations. 

The computer processed some 50 000 measurements and provided point-to-point, 35-mm-film 
plotting as well as printout of over 15 000 values of refractivity. The next section of this report pre- 
sents the surface-refractivity measurement results. 

SURFACE-REFRACTIVITY MEASUREMENT RESULTS 
Monthly Variations 

Table 2 summarizes the maximum, minimum, and monthly averages of refractivity as obtained 
from the 1 5 Apollo sites. The Antigua Island (ANG) station did not furnish pressure data; therefore, a 
nominal atmospheric pressure of 76.2 cm of mercury was used to obtain the ANG values in Table 2. 
This value of pressure corresponded to the average at Carnarvon, Australia, which is approximately 
the same height above sea level as Antigua (see Table 1). The maximum absolute variation in pressure 
at Carnarvon throughout the test was 1.5 cm of mercury. It is therefore estimated that the ANG 
values are correct to at least ±5 N-units. Also, it was noted (after the fact) that the GDS August and 
September 1967 pressure data were in error. Therefore, the values for these 2 months in Table 2 
were computed with P = 910 mb, based on 1968 GDS data; the results for these 2 months were not 
plotted. 

Figure 1 shows the measured total atmospheric average pressure during mid-May 1968. The 
expected correlation between station altitude and pressure (Figure 1) is quite apparent. In this test, 
the pressure values were seen to be primarily a function of station altitude, with maximum variation 
of only a few percent of total pressure throughout the test. 

The computer plots of refractivity measurements versus time from all sites (except ANG) are 
presented in Appendix A. These results clearly show the degree of variability in surface refractivity 
one can expect throughout any given month. The 1 -sigma spread is seen to be typically ±10 to 
±15 iV-units over a given month, depending on geographical location. The value of SD is a measure of 
the variability throughout a given month for a particular station. It is not an indication of measure- 
ment accuracy. The maximum and minimum values of refractivity reflect temperature extremes (low 
temperatures resulting in high refractivity and vice versa). The plots indicate the diurnal variations, 
which again are primarily temperature dependent, with lower temperatures generally being associated 
with local night. 
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Table 2-Monthly values of refractivity in iV-units. 
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Figure 1— Example of measured atmospheric pressure distribution. The pressures are 

averages for May 25 to 26, 1968. 


Monthly Averages Compared With Apollo Values 

Table 3 lists the monthly average refractivity estimates used in the Apollo program. Table 4 pre- 
sents the differences between the measured monthly averages and the corresponding Apollo estimates. 

In general, the agreement is good (within 13 percent or better for all cases), considering the fact that 
the measured monthly averages ranged from 266 to 385 in 1968, with individual values ranging from 
231 to 434. Where the Apollo refractivity averages are consistently above or below the observed 
monthly averages (for example, BDA or CRO), some adjustment of the Apollo values may be warranted. 

Table 5 compares the August 1967 to August 1968 and September 1967 to September 1968 
monthly average refractivity. The year-to-year variation in the monthly average refractivity in all cases 
is less than 6 percent and in several instances is less than 0.5 percent. 

CONCLUDING REMARKS 

The monthly average refractivity values used in the Apollo program are in all cases within 1 3 per- 
cent of the measured values. 

At any given site, humidity and temperature variations were seen to be the primary cause of 
changes in refractivity. Changes in air pressure at a given site (in the absence of temperature and humid- 
ity variations) were noted to result in a maximum change in refractivity of approximately ±5 A-units 
per year. 

This test has shown given site time variations in surface refractivity, a scaling quantity in most 
mathematical models for the troposphere. A necessary continuing study is the comparison of mathe- 
matical models of the troposphere with measured profiles (e.g., Reference 13). An efficient and accu- 
rate algorithm for computing satellite tracking data tropospheric corrections is given in Reference 14. 
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Table 3 -Monthly average refractivity estimates used in the Apollo program. 


Station 






Month 






Jan. 

Feb. 

Mar. 

Apr. 

May 

H 




El 


m 

ACN 

361 

364 

366 

363 

361 

356 

353 

353 

354 

356 

355 

357 

ANG 

360 

359 

359 

361 

372 

375 

378 

379 

376 

377 

372 

366 

BDA 

335 

333 

335 

343 

353 

372 

377 

378 

371 

359 

345 

338 

CRO 

338 

339 

341 

348 

362 

378 

385 

386 

381 

361 

349 

339 

CYI 

333 

336 

336 

335 

338 

346 

343 

352 

354 

350 

343 

336 

GBM 

343 

354 

352 

361 

370 

375 

386 

383 

384 

373 

364 

356 

GDS 

279 

281 

278 

279 

297 

277 

279 

279 

278 

277 

276 

279 

GWM 

372 

377 

378 

383 

381 

377 

373 

379 

381 

380 

375 

379 

GYM 

318 

322 

320 

321 

331 

351 

368 


372 

346 

328 

321 

HAW 

299 

298 

297 

299 

301 

305 

308 

310 

309 

309 

304 

302 

HSK 

306 

312 

309 

293 

300 

297 

296 

295 

298 

301 


301 

MAD 

295 

293 

295 

294 

298 

302 

299 



304 

299 

296 

MIL 

338 

339 

341 

348 

362 

378 

385 

386 

381 

361 

349 

339 

TAN 

310 

310 

313 

305 

298 

292 

291 

290 

291 

293 

303 


TEX 

328 

333 

340 

350 

362 

377 

395 

375 

372 

353 

336 

334 


Table 4-Observed minus Apollo values of N s * 


Station 

1967 

1968 

Aug. 

Sept. 

Feb. 

Mar. 

Apr. 

May 

June 

July 

Aug. 

Sept. 

ACN 

-17 

-16 

-15 

-15 

-14 

-14 

-13 

-16 

-19 

-17 

ANG 

— 

— 

— 

-8 

+ 1 

0 

0 

-3 

+1 

+2 

BDA 

-11 

-10 

-10 

-6 

-10 

-7 

-7 

-4 

-5 

-9 

CRO 

+9 

+9 

+ 19 

+ 13 

+ 19 

+ 11 

+ 14 

+13 

+ 15 

+21 

CYI 

-18 

-15 

-8 

-13 

-8 

-7 

-9 

-5 

-9 

-12 

GBM 

-4 

-3 

-20 

-14 

-6 

-1 

+9 

-2 

+2 

-2 

GDS 

+4 

+8 

+ 13 

+ 10 

+5 

-10 

+8 

+3 

+1 

-8 

GWM 

+7 

+4 

-9 

0 

-5 

+2 

+4 

+ 12 

+6 

+3 

GYM 

+ 1 1 

+ 1 

+7 

+2 

+6 

+9 

+ 17 

+ 14 

+ 14 

+6 

HAW 

+5 

+5 

+8 

+9 


+9 

+6 

+8 

+3 

+6 

HSK 

— | 


-24 

-14 

-2 

-12 

-17 

-17 

-17 

-20 

MAD 

-17 

-20 

0 

-3 

+ 1 

-9 

-15 

-14 

-34 

-13 

MIL 

-7 

-9 

-24 

-15 

-6 

-10 

-8 

-12 

-14 

-16 

TAN 

-3 

0 

-2 

-3 

-3 

-3 

-3 

+3 

+1 

-1 

TEX 

+ 1 

0 

+5 

+ 15 

+25 

+ 19 

+8 

-10 

+8 

+4 


*The Apollo values are NASA-MSC best estimates of monthly averages at each site. The same values are used each year. The range 
of the measured 1968 monthly average AN S is 1 19, the maximum being 385 (TEX, June 1968) and the minimum being 266 (MAD, 
August 1 968). These values are best estimates based on past meteorological surveys such as compiled in Reference 10. 


9 



























Table 5— Year-to-year variation in the August and September averages of N s for 1967 and 1968. 


Station 

August 

1967 

August 

1968 

lAVl 

September 

1967 

September 

1968 

lAA/l 

ACN 

336 

334 

2 

338 

337 

1 

ANG 

no data 

380 

— 

no data 

378 

- 

BDA 

367 

373 

6 

361 

362 

1 

CRO 

331 

338 

7 

334 

346 

12 

CYI 

334 

343 

9 

339 

342 

3 

GBM 

379 

385 

6 

381 

382 

1 

GDS 

283 

280 

3 

286 

270 

16 

GWM 

386 

385 

1 

385 

384 

1 

GYM 

381 

384 

3 

373 

378 

5 

HAW 

315 

313 

2 

314 

315 

1 

HSK 

no data 

278 


no data 

278 

- 

MAD 

283 

266 

; 7 

286 

293 

7 

MIL 

379 

372 

7 

372 

365 

7 

TAN 

287 

291 

4 

291 

290 

1 

TEX 

376 

383 

7 

372 

376 

4 
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Appendix A 


Computer Plots of Refractivity Versus Time 
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